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We provide a family of Non-Gaussian pure entangled states in a bipartite system as the eigen 
states of a quadratic Hamiltonian composed of the Einstein- Podolsky-and-Rosen-like operators. The 
ground state of the Hamiltonian corresponds to the two-mode-squeezed vacuum state, and belongs 
to the Gaussian states. In contrast, all of the excited states are Non-Gaussian states. A separable 
inequality maximally violated by the eigen states is derived. 
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The number states are fundamental states to describe 
quantum nature of the world whereas the coherent states 
are more accessible Gaussian states in quantum optics 
and of relevance to the Gaussian description of the na- 
ture. The single-photon state and more-than-one-photon 
number states are typical example of Non-Gaussian (NG) 
states and of essential for the quantum information pro- 
cessing (QIP) [H, 0]- In the combined systems, the en- 
tanglement property of the states are one of the central 
objective on the modern physics as well as on QIP [1, 
For the class of Gaussian states, the characterization of 
entanglement has been developed [3, M, 0, 0, S] , and there 
are several approaches to take into account the NG states 
i, i, [3 EH, d • However, it seems that the number of 
the examples of the NG entangled states is rather lim- 
ited and the testing ground for entanglement theories 
concerning the NG states is lacking. 

In this Letter, we introduce a family of NG pure en- 
tangled states as the eigen states of a Hamiltonian that 
describes two harmonic oscillators. The states are con- 
sidered to be a typical example of the NG states similar 
to the case of the number states in the one mode system, 
and serve as a useful tool to analyze the entanglement 
in infinite dimensional systems beyond the Gaussian ap- 
proach. 

In the seminal paper by Einstein, Podolsky, and Rosen 
(EPR) [3] , the nonlocal property of quantum continuous 
variable system was demonstrated by introducing simul- 
taneous eigenstates of the EPR operators Xj^ — xb and 
Pa + pB , those are uncommutable on each of the sub- 
systems — [xbtPb] = h but commutable on the 
total system, [xa — xbtPa+Pb] = 0. The essential point 
for nonlocality is that the total uncertainty is smaller 
than the sum of the subsystem's uncertainty, and such 
property can be seen for the EPR-like operators 



XA /T. 



PA , /T, 
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Associated with this pair of the EPR-like operators, 
Duan et al., provided the separable criterion [7|: Any 
separable state must satisfy 

m)) > r'+c 



Here the total noise of the EPR-like operators is defined 

by 



(3) 



and AO := O — (O) represents the deviation of the oper- 
ator. In what follows we assume (xa) — (xb) — (Pa) — 
(Pb) = and < ^ < 1. 
From the commutation relation (i_ (^)] = 

— the lower bound of (57) is shown to be 
|[a'+(^),tt'_(^)]| = This minimum of f2 is achieved 

by the two- mode-squeezed vacuum (TMSV): 
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where the number states are given by \n)A ■= 
(at)"|0)A/V^ and \n)B := (&^)"|0)B/Vn! associated 
with the annihilation operators of mode A and B, a :— 
{xA + ipA)l^2 and h :— (xb + ipB)/V^, the vacuum 
states defined by a\0)A — and 5|0)b = 0, and the 
commutation relations: [a, a^] = 1, [fe, 5^] = 1, and 
[a, 6] = [a, 6t] = 0. 

Using these ladder operators and their commutation 
relation we can write 



Let us define the new annihilation operator by 

a-0 



(5) 



This implies \A^, = 1 and f7(^) takes the form of the 
Hamiltonian of a single harmonic oscillator. 



m = (r'-e)(2ip? + i 



(2) 
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From this expression it is clear that (f2(^)) is no smaller 
than ((^~^ — ^) since A^ A > 0. The state that achieves 
the minimum of is also derived from the condi- 

tion A^\0,0;S^) = 0. From Eqs dU and (O we can ver- 
ify A^\'ip^)AB = 0, and thus denote |0, 0;^) := |V'5)as- 



2 



The eigen states of the number operator Na '■= IS 

generated by {AI)'^^\0,0;^) =: \Na,0;O and the 
eigen value Na is a non-negative integer. The eigen state 
\Na,0;O reduces to \Na)a\0)b in the hmit ^ -> 0. This 
suggests that there is another harmonic oscillator whose 
eigen-state reduces to \0)a\Nb)b in ^ ^ 0. We can find 
another pair of the EPR-like operators, i^±{^^^), which 
commute with the other ones as [w±(C)i^±(C "'^)] = 0- 
Associated with a)±(if~^) we define another annihilation 
operator, 



6 -Cat 
7^ 



(7) 



Then, we can write the sum of the total variance of the 
two pair of the EPR-like operators as 



1 .(8) 



Noting that B^|0,0;0 = and [A^^B^] = [A^,^] = 0, 
we obtain all the eigen states of + Cl{£^^^): 



\Na,Nb;0 = 



|0,0;0. 



(9) 



We call this state the entangled number state (ENS) 
or two-mode-squeezed number state 0, While the 

vacuum state |0, 0;^) is a Gaussian state, the excited 
states are NG states. This can be seen by combin- 
ing the orthonormal relation {N'j^, Ng;£^\NA, Nb;0 — 
^NA,N'^SNB,N'g with the fact that a state orthogonal to a 
Gaussian state is a NG state. 

Note that not only the TMSV but also the states 
\0,Nb;£,) with arbitrary excitation of Nb achieve the 
minimum of il(C). This implies that the minimum is 
also achieved by the infinite sequence of the NG states 
{|0, Nb; \-^B > 1} and these states are entangled. For- 
mally, the eigen space of Na specified by Na = is 
spanned by {\0, Nb; £,)\^b > 0}, and any state belongs 
to this eigen space achieves the minimum. From the sep- 
arable criterion of Eq. Q we can obtain a somewhat 
strong statement: any state spanned by {|-/Va, Nb; 0\0 < 
Na < CV(1 - ^^),Nb > 0} is entangled. 

In order to see that the ENSs are entangled state, 
we will write \Na, Nb;0) in the photon number basis 
{\n)A\m)B}- In what follows we assume Na > Nb for a 
simplicity. Using the relation e"^" ^ 6e^° ^ = b + ^a^ 
we have (ip^^e«'^'''' = (1 - e)~^^/^[a^{l - e) - 
^b]^"^. By acting this operator on the vacuum, we ob- 
tain \Na,0;O = (i|)^-v/T^e«"'^>)^|0>B/ViVlT- 
(1 _ C2)(iVA+i)/2g5at6t|^^^|Q^_ Similarly, from the rela- 
tion e-«'^^''^ae«'^*''^ ^a + ^b^ we have \Na,Nb;0 = (1- 



Then we find the Schmidt decomposed form of the ENS: 
\Na,Nb;0 



^^-y^\,aw ^ NbKi - e)H-ir^-'' 



k=0 



k\{NB ~ k)\ 



k\ 



Na\ 



-\NA-NB + k)A\k)i 



Nb\ \ (NA-NB + ky: 

OQ 

Y,Cm{NA,NB,e)\NA-NB+m)A\m)B (10) 



m=0 



where the Schmidt coefficient Cm is given by [14 1 
C™(iV^,7VB,C') 



m.'m{m,NB } / 
fc=0 \ 



^NA\NBKNA-NB+my.m\ 



fc!(m - k)\{NA -Nb + ky.{NB - fc)! 



(11) 



Physically, {|C„ip} means the photon number distribu- 
tion of the system B (or A) . Actually, we have the mean 
and variance of the number operator ni, := Wb, 



m 



m\Cr, 



{rib 



m=0 



NB+ejNA + i) 
1-e 



A^m := ^(m-m)2|C™|^ 



m=0 

(hi) - {hb 



2, ,2_ ^HNa + Nb + 2NaNb + 1) 



(12) 



Since the distribution has non-zero variance A^m > 
whenever ^ > 0, any of Cm is less than unity and thus 
the ENS cannot be a product state. Therefore, all of the 
ENSs are entangled states. 

The degree of entanglement for the pure states 
are characterized by the von-Neumann entropy of its 
marginal state [l^. The entanglement corresponds to 
the Shannon entropy for the distribution of the squared 
Schmidt coefficients: 



Ei\NA,NB;0) 



J2 |C™|2l0g2|C„ 



(13) 



It reduces to the entanglement of the TMSV when Na = 
Nb = [51]. 

Figure [T] displays E{\Na,Nb;0) for ^ = 0.7 and 
{{NA,NByO < Nb < Na < 10}. In order to grasp 
the entangled property of the ENS, it might be instruc- 
tive to consider the distribution {|CmP} for one- mode 
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FIG. 1: The entanglement for the entangled number states 
\Na,Nb;£,) with ^ = 0.7 and < iV^ < A*'s < 10. The state 
is specified by {Na,Nb)- 

excitations, e.g., we set Nb = 0. In these cases, {|Cmp} 
corresponds to the Negative Binomial (NB) distribution 

|C„(A^^,0,p)p = (1 -p)i+^V"i%±^, (14) 

and the reduced state of |Afy!i,0;^) is a NB state [l5| . 
As is shown in FIG. [2l the NB distribution has a single- 
peaked waveform whose mean is and variance is 
^\i-p)^^ . It is also known that the NG distribution ap- 
proaches to the Poisson distribution P(fc; A) = \^e~^ /k\ 
with A = [Na + 1)(1 —p)/p for larger Na- The excita- 
tion of Nb modulates the normal waveform of the NB 
distribution and adds the nodes as exemplified in FIG. 
[5J The parameter ^ determines the typical scale of the 
distribution [see also Eq. fT^ ]. and we can see that the 
distribution becomes broader as ^ — > 1. From these ob- 
servations on the distribution of |CmP and the fact that 
the entropy becomes higher if its distribution becomes 
more uniform, it is likely that the ENS has the following 
property of entanglement: (i) The entanglement of ENS 
becomes larger when Nb becomes larger with fixed Na] 
(ii) The entanglement of ENS becomes larger when Na 
becomes larger with fixed Nb] and (iii) The entangle- 
ment of ENS becomes larger when ^ becomes larger with 
fixed Na and Nb- Similar conjectures have also been 
presented in 

The ENS is an eigen states of the number operator 
A^A^, and the variance (A^(A^yl^)) is zero. Hence, it is 

likely that the state who has a small variance of A\A^ 
is also entangled. To find out such a relation, let us 
consider the partial transposition (PT) [4] of the vari- 
ance, {(A^^A{)'^) PT - {MA^)%j,. Calculating the PT 
by using the relation ^ {a^^a^a'^''''a'^h'^%'Ph'^W) pt = 



FIG. 2: The distributions of the squared Schmidt coefiicient 
|CmP of the entangled number states \Na,Nb\0 with ^ = 
0.7. The state is specified by {Na,Nb). For Nb = 0, the 
coefficient |CmP is a normal wave packet whose mean is {Na + 
1)^/(1-^^) and variance is (iVA + l)C^/(l-C^)^. An excitation 
of the Nb mode induces a node on the distribution as is shown 
in the case of Na = 120 and A'^b = 0, 1, 2, 3, 4, where |C,„ p + 
0.02A''b is displayed instead of |Cmp. 

(at'=a'at"a"(6to6P6t95'-)t) ^ {a'f''a^a^"'a'^b^''b%^Pb°) we 
have the inequality 

{i-er {{iAiA^r)PT-{AiA^fpT) 
= {A'[{a^-^b^){a-m)+e>e, (15) 

which holds for any valid density operator. The PT of 
the density operator is also a valid density operator for 
separable states. From this fact with Eqs. ([6]) and (fT5|) 
we have a separable criterion: Any separable state must 
satisfy 

This inequality is maximally violated by the ENSs, which 
yield {A'^[Cl{^)]) = 0, and thus ENSs are negative PT 
entangled states. For a given Na, any state spanned 
by {\Na, Nb]0\^b > 0} also maximally violates the 
inequality (fTB|) . Interestingly, the two separable criteria 
of Eqs. ([2]) and (flB)) are concerned with quite different 
aspects of the same operator of Eq. ([3]). While 

Criterion ^ is of relevant to the Gaussian states our 
Criterion is, at least, capable of detecting the eigen 
states of Eq. ([5]), which form a family of the complete 
orthonormal set specified by the parameter ^. 

The ENS is an eigen state of the quadratic Hamilto- 
nian of Eq.®, and thus the lower exited states as well as 
TMSV will appear in such physical systems. By not- 
ing that the vacuum \0)a\0)b and TMSV |0,0;0 are 
connected with the unitary operator of the parametric 
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amplification U = e^"^''^ as |0,0;0 = U\0)a\0)b 
with r = tanh^^^ and from the relation j4^|0,0;^) — 
Ua^fWU^A^B = Ua^O)\0), the ENS will appear as 
the output of the parametric amplification when the in- 
put state is a number state Q, The operator 
A' (x a'^ + implies the process where the addition of a 
photon on the mode A and the subtraction of a photon in 
the mode B are done coherently. Hence ENS might also 
be realized by performing such process on the TMSV. 

Given the "number-state" basis we may define the "co- 
herent state" or "squeezed states" . We define the "coher- 
ent state" as the simultaneous eigen state of the annihila- 
tion operators: A^\a,l3;^) — aja,/?;^) and I3^\a,P;^) — 
/3|a,/3;f) with complex amplitudes a and f3. It can be 
represented in the basis of the ENSs as 

= D^Ja)^Bj/3)|0,0;O, (17) 

where we used the relation e""*~"*° = g-l"l^/2gaa*g-a*a 
and introduced the displacement operator 

DA,{a):=e''4-'-'^^. (18) 

Since the displacement operator is written by the product 
of the local displacement operators the "coherent states" 
can be converted with each other by the local unitary 
operation and we have 

10,0.0. 

Hence, the "coherent states" are equivalent to the 
TMSV under the local unitary transformation, and 
the amount of entanglement is equal to that of the 
TMSV. The displaced TMSV states are also called as 
the two-mode-squeezed coherent states. Note that we 
have the product of the coherent states \a,(3;£,) — > 
\ct-)A\P)B in the limit ^ ^ 0. Note also that we 
can use the expansion by the over complete rela- 
tion with the equally entangled Gaussian pure states 
/ |a, f3; ^(q;, /?; £,\cP acP f3 / tt'^ = 1^® Is whereas we have 
the expansion by the complete orthonormal basis of the 
ENSs each of which has typically different amount of en- 
tanglement T,N^,Ns 1^^' Nb;6{Na, NB;^\ = tA(g)tB. 
These complementary roles played by the number states 
and coherent states are also thought to be useful in the 
theory of entanglement. 

In conclusion we have introduced the entangled num- 
ber states as the eigen states of a quadratic Hamilto- 
nian. The ground state corresponds to the two-mode- 
squeezed vacuum and belongs to the Gaussian states. 
The excited states are Non-Gaussian entangled states. 
The Schmidt coefficients and the entanglement were cal- 
culated. It was conjectured that the higher number state 



has larger amount of the entanglement. We have also 
derived a separable criterion from the partial transposi- 
tion of the variance of the quadratic form. This sepa- 
rable inequality is maximally violated by the entangled 
number states and these states are shown to have neg- 
ative partial transposition. With their mathematically 
tractability and well-defined properties due to the eigen 
states of the harmonic oscillators, the entangled number 
states could form a useful testing ground for entangle- 
ment theories concerning Non-Gaussian states, and the 
characterization of the entangled number states will be 
definitive starting point for understanding the class of 
Non-Gaussian entangled states. 
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